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SUMMARY 


The inverse problem for jet acoustics, or the determination of noise sources from far-held 
pressure information, is proposed as a tool for understanding the generation of noise by 
turbulence and for the improved prediction of jet noise. An idealized version of the problem 
is investigated first to establish the extent to which information about, the noise sources 
may be determined from far-held pressure data and to determine how a well-posed inverse 
problem may be set up. Then a version of the industry-standard MGB code is used to 
predict a jet noise source spectrum from experimental noise data. 
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1. INTRODUCTION 


The purpose of this paper is to describe an initial examination into the use of the inverse 
problem in acoustics to understand properties' of turbulence-generated noise. That is, it is 
desired to use experimental determinations of the far-field pressure to generate information 
about the generation of noise by turbulent flows. Naturally, such an inverse problem must 
be carefully set up in order to give a mathematically well-posed problem. If this may be 
done successfully, however, there is a potential for eliciting important information about the 
still poorly understood process of the generation of noise by turbulence. 

There has recently been a number of investigations into the use of sound waves as a 
diagnostic tool for fluid flows [1, 2, 3]. This technique involves passing sound waves through 
a flow and reconstructing the flow based on the changes in amplitude and phase in the 
sound waves produced by velocity variations in the flow. What is proposed here is somewhat 
different, both in technique and emphasis, because here we wish to learn about the nature of 
the sound sources rather than the moving fluid through which the sound propagates and the 
inversion is based on the sound waves produced by the flow itself, rather than sound waves 
applied externally. 

The analysis of the inverse problem in jet acoustics discussed in this paper is based on 
the Lighthill acoustic analogy [4] and its refinements [5, 6]. We begin with an examination 
of the inverse problem in its simplest context, where the basic Lighthill acoustic analogy is 
used to relate far-field pressure information to the noise sources. This serves to establish 
how a well-posed inverse problem should be set up and what information we can and cannot 
expect to get from solving the inverse problem. 

Having established the basic features of the inverse problem, we apply the MGB code 
of Khavaran [7] as modified in previous work [9] to permit arbitrary noise spectra to be 
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specified. Incorporating various corrections for convective and other effects, the MGB code 
provides a model for those aspects of the noise source that cannot be inferred from the 
far-field pressure data. We begin by attempting to determine the frequency distribution 
of the noise source from experimental data, with the remainder of the source information 
specified by the MGB code. Experimental data presented by Seiner et al. [10] is used in 
conjunction with aerodynamic data produced in earlier work [11] as inputs for the inverse 
problem solution procedure. 


2. AN IDEALIZED INVERSE PROBLEM 


In this section, the basic Lighthill acoustic analogy is used to relate noise sources to the 
far-field pressure spectrum in a relatively simple way. Idealized representations of the jet 
noise problem are used to determine how a well-posed inverse problem may be set up and to 
determine what information can be derived from the inverse problem and what information 
cannot. The experience gained through the examination of this idealized inverse problem 
will be applied to the formulation and numerical solution of the inverse problem for a jet 
flow in the following section. 

Lighthill advanced an acoustic analogy [4] in 1952 for the analysis of flow-generated 
noise. He recast the Navier-Stokes equations as a wave equation for the acoustic pressure 
disturbance p that includes source terms involving the nonlinear convective terms, 


1 d 2 p 
cl dt 2 


- V 2 p = 


dxidxj 


(1) 


The quadrapole source strength T t j is typically approximated [7] by only the term quadratic 
in the velocities, the viscous stress and other terms being deemed neglegible in the linear far 
field: 


Tij Po^i^j . 


( 2 ) 
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Here, the fluid density has been approximated by its ambient value p„ and the Vi are the 
velocity components. The solution to this equation may be expressed using the frec-space 
Green’s function for the three-dimensional wave equation as 


p(x, t) 


d 2 f, Tij{ y,t- [x-y|/c 0 ) 
dxidxj J 47r|x — y| 


( 3 ) 


At distances large compared to the linear extent of the source’s support, this expression 
may be approximated by 

P(M) = 4^x1 / dy Ttj( - y,t ~ |X ~ y|/Co) (4) 


The quantity most commonly measured in acoustics experiments is the far-field mean-square 
pressure; it may be represented now as 


(p(x, t)p(x, t + r))= / d y 7 o«(y> T ) ( 5 ) 

where 

Wy.O = |^/df<iyzi,). (6) 

Here, Tjj is evaluated at position y and time £, T' kl is evaluated at y + ^ and t + r; ^ is the 
spatial and r is the temporal separation of the two source strengths T tJ . The statistics of 
the source strengths are assumed to be time stationary. 

The two spatial integrations in the term on the right-hand side of (5), over y and £, have 
the effect of hiding information about the sources from the far-field pressure and so there is 
no hope of getting complete information about the noise sources without instituting a special 
experimental program specifically set up to probe the nature of the sources, perhaps along 
the lines of References [1, 2, 3]. While this is an intriguing possibility that is discussed further 
in the conclusion, the purpose of the present investigation is to determine what information 
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about the noise sources may be derived based on current standard experiments and existing 
experimental data. 

Typically, experimental data is expressed in terms of the Fourier transform of (5) with 
respect to r and so, in practice, in the following, we will be concerned with attempting to 
determine the noise spectrum Iijki(ty- Furthermore, because experimental data is generally 
available only on some arc in the far field, parametrized by some angle, say 9, the data 
available to us for solving the inverse problem is a function of the two variables H and 9 
and we cannot hope to fully determine the dependence of f dylijki on x. Rather, the inverse 
problem as posed here has the potential for charaterizing one temporal and one spatial degree 
of freedom of the noise source. The remainder of the noise source function must be specified 
a priori. 

Additional information about the source term on the right-hand side of (5) may be 
deduced using the fact that (T^T^) is essentially the fourth-order correlation of a turbulent 
velocity field. Thus, the characteristic length and time scales for the separation distance £ 
and time r are given by the correlation lengths and times of the turbulence. On the other 
hand, the variable y in I^i gives the position of the source in the flow field, and gives the 
dependence of the turbulence correlation on the aerodynamic data through such quantities 
as the kinetic energy and the dissipation. The fact that both y and £ are integrated over 
in determining the pressure in (5) means that, without additional experimental data, only 
limited information about the spatial dependence of the turbulence correlations and the 
aerodynamic data may be inferred when the assumptions leading up to (5) are satisfied. 
Consequently, we focus on determining information about the temporal dependence of the 
correlation (the frequency spectrum) and the directivity (the tensor nature of Iijki)- 

3. THE INVERSE PROBLEM VIA THE MGB CODE 
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The MGB approach to noise prediction [8] provides a convenient basis for an initial 
analysis of the inverse problem in a real jet problem. In the present work, use is made of the 
MGB code developed by Khavaran [7], which employs the Lilley formulation of the Lighthiil 
acoustic analogy [5] and convective corrections as proposed by Ffowcs-Williams [6]. This code 
has been modified to permit alternative noise source spectra to be tested, as described in an 
earlier paper [9]. Here, this modification will be used to implement an iterative procedure, 
based on the MGB code, for the determination of a noise source spectrum that yields a 
given, experimentally determined, far-held noise spectrum. 

The MGB code, as modified in [9], approximates the noise source spectrum in the Lighthiil 
analogy (1) by the expression 

hjkiity — hjki 7 ( 0 ), ( 7 ) 

where 

904R7T 2 r°° o r+oo 

7(0) = — p 2 fi 4 / dk [kQ{k)f / [r{k,r)} 2 e inT dr (8) 

15 J 0 J-OO 

and 

lijkl — ^ + fiil$kj) + 25ij(ifcz] . (9) 

These expressions result from the standard assumption in the MGB analysis that the fourth- 
order correlations in Iijki may be written in terms of second-order correlations and the writing 
of those second-order correlations Sij in a form appropriate for statistically homogeneous and 
isotropic turbulence [12]: 


Sij&r) = Jdk e^Piji k) Q(k) r(k,r). (10) 

Here, the tensor Pi 3 { k) = 6^ — kikj/k 2 ensures that 5 \ 3 is isotropic and solenoidal, Q(k ) is 
computed from the energy spectrum 


Q(k) = 


1 

4irk 2 


m 


(ii) 
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and r(k , r) describes the temporal correlation at different wavenumbers; it is normalized so 
that r(k, 0) = 1. 

The discussion of the previous section indicates that we can expect to obtain via the 
inverse problem information about the frequency spectrum /(ft) and the directivity Iijki] 
the various assumptions about the turbulence statistics listed above serve to eliminate the 
remaining uncertainties in the source strengths. Once this information is obtained, some 
properties of the energy spectrum Q(k ) and time correlation r(k,r ) may be inferred, but 
clearly they cannot be completely determined. Naturally, any conclusions drawn are predi- 
cated on the correctness of the assumptions leading up to this representation of the source 
strength. 

We concentrate first on obtaining the frequency spectrum /(ft) from experimental data; 
determination of the directivity Tij k i will be dealt with later. Experience in obtaining curve 
fits for the spectra discussed in Reference [9] suggests that an appropriate representation of 
the frequency spectrum is 

/(ft) = exp [-E^.jOift 1 ] - (12) 

The coefficients a t were determined using a Newton-type method, with the MGB prediction 
code performing the function evaluation. The difference between the target far-held pressure 
and that predicted with the MGB code and given values of the coefficients was driven 
to zero by determining the “design sensitivities” with respect to the “design variables” Oj 
and solving for the increments to the a,. The necessary derivatives were computed using 
finite differences; that is, each of the a, was perturbed in sequence to give the columns in 
the sensitivity matrix. This simple scheme proved to require quite good initial guesses for 
the coefficients in order to converge and so a continuation scheme was implemented in which 
a number of intermediate target far-held pressure distributions were computed, gradually 
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moving towards the experimental target. 


4. RESULTS 

The MGB-base procedure described in the previous section was applied to the case of the 
Seiner jet [10] at 1550°F. The acoustic data employed in the inversion is that of Reference 
[10] and the aerodynamic data is from Reference [11]. The relevant plot from the latter 
reference is reproduced here as Figure 1. The ISAAC results with the k — e model, but with- 
out a compressibility correction, seemed to give the best agreement with the experimental 
aerodynamic data and are used in the present analysis. 

For an initial test of the concept, data collected from a microphone at an angle of 93.3° 
from the jet inlet was used to construct a three-parameter approximation to the source 
spectrum 1(0,). (Of the coefficients in (12), a 0 , Oi and 02 were nonzero, the others were 
zero.) The results are shown in Figure 2 where it is seen that a good fit is possible and that 
the present results are significantly better than those of the standard MGB model for the 
source spectrum. 

The source spectrum determined through the inverse problem is shown in Figure 3, along 
with the MGB frequency spectrum for comparison. The present spectrum is dramatically 
different: it has its peak at a much lower frequency and seems to have a fatter tail. The 
failure of the of the present spectrum to tend to zero as the frequency becomes small — as 
it must on physical grounds — is due to the absence of low-frequency experimental data in 
the target pressure distribution: if the pressure data had dropped to zero, then so would the 
result of the inverse computation. 

The most important conclusion to be deriv ed from this result is that the frequency 
spectrum of the noise sources seems to have a much longer time scale than that employed by 
the standard MGB code. As the MGB time scale is inherited from the second-order velocity 
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correlations through the analysis outlined in the previous section, this would seem to imply 
that the fourth-order correlations of the turbulent velocity field have a longer correlation 
time than the second-order correlations. Whether this is in fact the case or only seems to 
be in the present context as a result of some acoustic filtering effect akin to that described 
by Wilson et al. [13] is currently under investigation. 
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LIST OF FIGURES 


Figure 1: 1550° F jet case, no compressibility correction. Experiment, solid line; ISAAC k — 
e model, dash-dot-dot line; ISAAC GS ASM, long-dash line; PAB3D k — e model, 
dotted line; PAB3D GS ASM, dash-dot line; PAB3D Girimaji ASM, dashed line. From 
Reference [11]. 

Figure 2: SPL at a microphone angle of 93.3° for the 1550°F jet. Experiment, solid line. 

Khavaran MGB code, dashed line. Present compution, using source spectrum from 
inverse problem, dotted line. 

Figure 3: /(U)// max for the the noise spectrum found by solving the inverse problem (solid line) 
and for the MGB noise spectrum (dashed line). 
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Figure 1: 1550° F jet case, no compressibility correction. Experiment, solid line; ISAAC 
k — t model, dash-dot-dot line; ISAAC GS ASM, long-dash line; PAB3D k — e model, dotted 
line; PAB3D GS ASM, dash-dot line; PAB3D Girimaji ASM, dashed line. From Reference 
[ 11 ]. 
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Figure 2: SPL at a microphone angle of 93.3° for the 155CPF jet. Experiment, solid line. 
Khavaran MGB code, dashed line. Present compution, using source spectrum from inverse 
problem, dotted line. 
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